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We simulate competitive two-component growth on a one dimensional substrate of L sites. One 
component is a Poisson-type deposition that generates Kardar-Parisi-Zhang (KPZ) correlations. The 
other is random deposition (RD). We derive the universal scaling function of the interface width for 
this model and show that the RD admixture acts as a dilatation mechanism to the fundamental time 
and height scales, but leaves the KPZ correlations intact. This observation is generalized to other 
growth models. It is shown that the flat-substrate initial condition is responsible for the existence 
of an early non-scaling phase in the interface evolution. The length of this initial phase is a non- 
universal parameter, but its presence is universal. We introduce a method to measure the length of 
this initial non-scaling phase. In application to parallel and distributed computations, the important 
consequence of the derived scaling is the existence of the upper bound for the desynchronization 
in a conservative update algorithm for parallel discrete-event simulations. It is shown that such 
algorithms are generally scalable in a ring communication topology. 
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I. INTRODUCTION 

In (1+1) dimensions, the roughness of a surface that 
grows on a one dimensional substrate of length L can be 
expressed by the interface width wit) at time t 
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where hk{t) is the height of the column at site k and 
h{t) is the average height. The angular brackets denote 
the average over many configurations and the bar over a 
symbol denotes the average over L sites. The self-affined 
roughness of the interface manifests itself by the presence 
of the Family- Vicsek (FV) scaling 0: 
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where the scaling function f{y) describes two regimes of 
the width evolution: 
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The dynamic exponent z gives the evolution of the lat- 
eral correlation length ^(t) ^ t^^^. When £_{t) exceeds 
the system size L the width saturates. At saturation, 
for t 3> tx, the width scales as w ^ L", where a is 
the roughness exponent. The growth phase is the initial 
phase for t ^ ty before the cross-over time ^ 
to saturation. The growth phase is characterized by the 
single growth exponent f3 — a/z. The roughness, growth 
and dynamic exponents are universal. Their values de- 
pend only on the underlying mechanism that generates 
the correlations and scaling. 



A simple continuum model of nonequilibrium growth 
that leads to the scaling of a noise-driven interface is 
provided by the Kardar-Parisi-Zhang (KPZ) equation 0|. 
In the co-moving frame, the KPZ equation is 
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where h = h{x, t) is the height field (subscripts de- 
note partial derivatives). Coefficients ly and A give the 
strength of the linear damping and the coupling with 
nonlinear growth, respectively. The uncorrelated Gaus- 
sian noise Ci^it) has zero mean and covariance 



(C(x, t)C{x',t')) ^D5{x~ x')5{t - t'), 
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where D is the noise strength. A renormalization group 
analysis [2, |^ can provide a connection between the 
stochastic growth equation and scaling exponents. The 
KPZ universality class, governed by dynamics given by 
Eq. 10}, is characterized by a = 1/2 and fi = 1/3. A 
characteristic signature of the KPZ scaling is the expo- 
nent identity a -\- z = 2, valid in all dimensions. When 
A = in Eq. Q, the growth is governed by the linear 
Edwards- Wilkinson (EW) equation Q . The EW univer- 
sality class is characterized by a = 1/2 and /9 = 1/4. 
Using scaling arguments, it can be shown that in (1+1) 
dimension the EW exponent identity is 2a + 1 = z . 
When A = and v = Q va. Eq. 10}, the growth belongs 
to the RD universality class, characterized by (3 = 1/2 
and the lack of saturation. The RD interface is not self- 
affined. 

The theory behind kinetic rougheningand the origins 
of scale invariance are well understood 1^ IE IE S lE l^i 

but there are numerous instances of growth processes 
that neither follow one power law nor exhibit a clear-cut 
universality as it is expressed by the FV scaling. One 
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group of examples is the anomalous roughenin g in epi- 
taxial growth models 0, 0, 0| , fractures Il3j and 
in models with subdiffusive behavior or quenched disor- 
der 14]. These systems exhibit different dynamic scaling 
on local and global scales, characterized by different val- 
ues of roughness exponents. The super-rough dynamics 
of tumor growth is the first experimental observa- 
tion of anomalous scaling in (1-1-1) dimension. Another 
issue is the clear experimental observation of the KPZ 
universality and the role of quenched noise in the asymp- 
totic KPZ scaling [ElEEl- For one dimensional KPZ 
growth, by applying a weak noise canonical phase-space 
method, it has been shown recently that the KPZ dy- 
namic exponent is associated with the soliton dispersion 
law T?]. However, at saturation all KPZ correlations 
are exactly the same as would result from the linear EW 
equation 0, 0| . The fact that the EW equation is 
sort of "embedded" in the KPZ equation may give rise to 
ambiguous values of scaling exponents for growth mech- 
anisms (or models) that interpolate between the weak 
and the strong nonlinear coupling regimes. Even when 
the nonlinear coupling is strong the discrete models re- 
quire sufficiently large len gth and time scales to show 
clear KPZ scahng |ll IHIII IH Likewise, large 

scales are essential in simulation studies of roughening 
in the two-component growth models that combine one 
process governed by linear EW dynamics with another 
process governed by nonlinear KPZ dynamics |23|. Re- 
cently, Chame and Reis 2^ simulated in (1-1-1) dimen- 
sion a mixed growth where particles aggregated either by 
ballistic deposition (with probability p) or by random de- 
position with surface relaxation (with probability I ^ p). 
They show that for small p and sufficiently large L, the 
interface width has three well-defined evolution stages. 
The first stage, for early times, is the EW growth. The 
second stage, for intermediate times before saturation, is 
the KPZ growth. The third stage is the saturation. 

Two distinct growth phases were also observed in ex- 
periments with interfacial roughening in Hele-Shaw flows 
|26l |2^ |28| , in simulations of electrophoretic deposition 
of polymer chains |2^ [s^l and in numerical studies of 
one dimensional restricted solid-on-solid models with two 
growth components, each of which being a mechanism 
ruled by dynamics that belon gs t o a distinct universality 
class [IlllililllllilllllHllsl. The latter exam- 
ples suggest that the two-phase growth may be an effect 
of mixing the universalities. In support of this hypothe- 
sis comes a two-component model considered by d a Silva 
and Moreira [s^, where the Kim-Kosterlitz "wj depo- 
sition occurs with probability p and ballistic deposition 
occurs with probability (1 — p). Both of these growth 
mechanisms are governed by the KPZ dynamics, except 
that their corresponding continuum equations differ in 
the sign of the coupling A. In this case, large scale simula- 
tions in (1+1) dimension produce only one grow th p hase 
with the growth exponent /3 parametrized by p |39j . 

The purpose of roughening studies presented in this ar- 
ticle is to uncover the dynamics of desynchronization in a 



conservative parallel discrete-event simulations (PDES). 
The PDES are a technical tool to uncover the dynamics 
of information-driven complex stochastic systems. Their 
wide range of applications in contemporary sciences and 
technology 41] has made them an active area of research 
in recent years. Parallel and distributed simulation sys- 
tems constitute a complex system of their own, whose 
properties can be uncovered with the well-established 
tools of statistical physics. 

In PDES physical processes are mapped to logical 
processes (assigned to processors) that manage state- 
updates of assigned physical subsystems. The main chal- 
lenge arises because logical processes are not synchro- 
nized by a global clock. Consequently, to preserve causal- 
ity in PDES the algorithms should incorporate the so- 
called local causality constraint [i^. li^ ] . Depending on 
the way the local causality constraint is implemented, 
there are two broadly defined classes of u pda te proto- 
cols conservative al gori thms [i^. lil. l45l llq and 

optimistic algorithms |4^ l48l li^ . Important efficiency 
considerations of these algorithms involve the utilization 
of the parallel processing environment (as measured by a 
fraction of processors working simultaneously at a time) 
and the memory per processor required by state savings. 
The latter is closely related to the statistical spread in 
the processors' local times, i.e., to desynchronization. 
Recent applications of conservative PDES to modeling 
complex physics systems include ballistic deposition |50j ] , 
Ising spins 51] and dynamic phase transition |53.l53|. In 
an application to simulating Ising spins, an onset of self- 
organized critical behavior in optimistic PDES has been 
recently reported 0, [s^ . 

Since the introduction by Korniss et al. [5d| an idea of 
utilizing nonequilibrium surface growth methods 0j IS 
in evaluating the scalability of algorithms for conserva- 
tive PDES, there have been a number of advances using 
such efforts. The main concept behind this idea is the 
virtual time horizon (VTH) of the algorithm. The VTH 
is a (1-Fl) dimensional nonequilibrium surface. Its time 
evolution can be simulated by applying a deposition rule 
that is defined by a parallel-update protocol based on the 
algorithm. Several properties of the algorithm can be de- 
duced from analyzing its corres ponding simulated VTH. 
One of them is the utilization [53, ISg • Another one is 
the desynchronization of processors in the system as the 
PDES evolves in time. The width of the simulated VTH 
provides a measure of this desynchronization. It has been 
demonstrated that an asymptotic lack of synchroniza- 
tion 'SOf] can be avoided in new generations of algorithms 
(60b t^]' The focus of past studies was on the worst- 
case performance scenario when each parallel process 
consisted of only one computational object (or computa- 
tional volume). Despite being informative (e.g., provid- 
ing the evidence of a lower bound on the utilization) such 
studies are of little practical value since taking full advan- 
tage of parallelization implies many computational ob- 
jects per processor. Past simulation studies have demon- 
strated that when the load per processor is increased the 
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utilization increases sie nificantly IHil H l62| . and a 
mean-field-like method to estimate the utihzation in this 
general case has been developed Ell- The missing 
element thus far has been a detailed investigation of the 
dynamic scaling properties of the simulated VTH for the 
general case of many computational objects per proces- 
sor. This is the main theme of this article. 

We simulate three growth mechanisms for the VTH 
that correspond to three implementations of a conserva- 
tive update algorithm in a ring communication topology, 
where each processor communicates only with its imme- 
diate neighbors. The models are variations of Poisson- 
random depositions. Two of these rules are realizations 
of the worst-case scenarios and present a situation when 
the system of L processors can be mapped onto a closed 
chain of L spins. The third model corresponds to the gen- 
eral case when each processor carries N computational 
volumes and it combines the second model with random 
deposition (RD). In this work, we do not attempt to ob- 
tain exact values of scaling exponents that characterize 
these growth processes. Our primary interests are in the 
scaling functions for the interface width and in the uni- 
versal properties of the VTH interfaces. 

Simulation studies of the VTH interface in the worst- 
case scenario showed that for very lar ge L this interface 
belongs to the KPZ universality class jSg- However, for 
small L or for early times, before the KPZ growth is at- 
tained, the width does not scale. This suggests the strong 
sensitivity of the evolution to the initial condition. When 
the model is generalized to accommodate many computa- 
tional volumes per processor, the evolution of the VTH 
width changes. Now there are two distinct phases in 
the growth regime. The early phase evolves in the RD 
fashion and the later phase has signatures of the KPZ 
scaling. In this work we investigate the above issues in 
large-scale simulations. Unless stated otherwise, configu- 
rational averages were obtained over an ensemble of 800 
independent simulations. The VTH models and defini- 
tions are explained in Sec. ^ Simulations are initiated 
from a flat substrate and carried on up to 10^ time steps, 
well beyond cross-over times to the steady state for the 
considered substrate sizes. In Sec. Illll we analyze the in- 
terface evolution for random depositions at local surface 
minima (i.e., the worst-case scenario) and show that the 
initial lack of scaling is an artifact of the flat-substrate 
initial condition. To identify nonuniversal features in the 
evolution, in addition to Poisson-random depositions we 
also consider both Gaussian and uniform-random depo- 
sitions. In the steady-state time averages we omit the 
index t in the notation, e.g., {w^) denotes the saturated 
surface width. In Sec. llVl we perform the analysis of the 
interfaces generated by two simultaneously acting growth 
mechanisms, one of which being RD and the other gen- 
erating the KPZ correlations, and we obtain a universal 
scaling function for this type of VTH interfaces. Results 
obtained in this section show that the RD admixture 
elongates the principal height and time scales, leaving 
the KPZ correlations intact. In Sec. 0we derive a gen- 



eral relation between the VTH interface velocity and the 
utilization in conservative update processes. In Sec. lVl Al 
we generalize findings of Sec. lIVI to two-component mod- 
els that mix RD with a deposition that either classifies 
within the KPZ or within the EW universality class. In 
particular, we show that the RD admixture that hap- 
pens with probability (1 — p) gives rise to a p-dependent 
affine component in the scaling. Section IVIBI contains 
the discussion of finite-size effects observed in scaling 
of the VTH interfaces. An example of false scaling in 
Sec. IVI CI is provided to illustrate the importance of the 
relaxation from the flat-interface initial condition in the 
scaling considerations. Applications to scaling and scal- 
ability of conservative PDES algorithms are discussed in 
Sec. IVI Dl Conclusions are summarized in Sec. IVIII 



II. SIMULATION MODELS 

In simulations a system of L processors is represented 
as a set of equally spaced lattice points /c, /c = 1, 2, L. 
Each processor performs a number of operations and 
enters a communication phase to exchange informa- 
tion with its immediate neighbors. This communica- 
tion phase, called an update attempt, takes no time in 
our simulations. In this sense we simulate an ideal sys- 
tem of processors (the relation to PDES is discussed in 
Sec. IVI D|) . An update attempt is assigned an integer 
index t that has the meaning of a wall-clock time (in 
arbitrary units). 

The local virtual time hk{t) at the fcth processor site 
represents the cumulative local time of all operations on 
the fcth processor from the beginning at t = to time 
t. These local processor times are not synchronized by 
a global clock. The ring communication topology among 
processors is mapped onto a lattice arrangement with pe- 
riodic boundary conditions, ft.L+fe(t) = hk{t). The set of 
local virtual times hk{t) forms the VTH at t. The growth 
of the VTH is simulated by a deposition rule, where lo- 
cal height increments rjk (t) are sampled from the Poisson 
distribution of unit mean. The form of the deposition 
rule depends on the processor load, as explained below. 

A general principle that governs the conservative up- 
date protocol requires a processor to idle if at the up- 
date attempt t the local causality constraint may be vi- 
olated. This happens when at t the fcth processor does 
not receive the information from its neighboring proces- 
sor (or processors) if such information is required to pro- 
ceed in its computation. This corresponds to a situation 
when the local virtual time hk{t) of the fcth processor is 
ahead of either one of the local virtual times hk~i{t) or 
hk+i(t) of its left and right neighbors, respectively. In 
this unsuccessful update attempt the local virtual time 
hk{t) is not incremented, i.e., the fcth processor waits: 
hk{t -|- 1) = hk(t). In another case, for example, when 
at t the fcth processor does not need information from its 
neighbors it performs an update regardless of the relation 
between its local virtual time and the local virtual times 
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on neighboring processors. 

One example of computations that foUow the above 
model is a dynamic Monte Carlo simulation for Ising 
spins. In a parallel environment, a spin lattice is spatially 
distributed among L processors in such a way that each 
processor carries an equal load of one contiguous sublat- 
tice that consists of TV spin sites (i.e., each processor has 
a load of N volumes). Some of these N spin-lattice sites 
belong to border slices, i.e., at least one of their immedi- 
ate neighbors resides on the sublattice of a neighboring 
processor. Processors perform concurrent spin-flip oper- 
ations (i.e., increment their local virtual times) as long as 
a randomly selected spin-site is not a border site. If a bor- 
der spin-site is selected, to perform a state update that 
is consistent with and faithful to the underlying physi- 
cal spin dynamics, a processor needs to know the current 
spin-state of the corresponding border slice of its neigh- 
bor. If this information is not available at the t update 
attempt (because the neighbor's local time is behind), 
by the conservative update rule the processor waits until 
this information becomes available, i.e., until the neigh- 
bor's local virtual time catches up with or passes its own 
local virtual time. 

The least favorable parallelization is when each proces- 
sor carries one computational volume, iV = 1. Computa- 
tionally, this system can be identified with a closed spin 
chain where each processor carries one spin-site. At each 
update attempt each processor must compare its local 
virtual time with the local times on both of its neigh- 
bors. 

The second least favorable arrangement is when pro- 
cessors have a computational volume N = 2. As be- 
fore, the system can be mapped onto a closed spin chain 
where each processor carries two spin-sites, each of which 
is a border site. At each update attempt every processor 
must compare its local time with the local time of one of 
its neighbors. 

In general, when iV > 3, at update attempt t, the 
comparison of the local virtual times between neighbors 
is required only if the randomly selected volume site is 
from a border slice. In all cases, we start the simulation 
from a flat substrate at i = 0, hk{Q) = 0. 

At every successful update attempt, the simulated lo- 
cal virtual time at the kth site is incremented for the 
next update attempt: hk(t + 1) = hk(t) + rik{t), where 
Vkit) = — In(rfet), and r^t G (0; 1] is a uniform random 
deviate. The three cases described above are realized in 
simulations by the following three deposition rules. 
Rule 1 [N = 1): The update attempt at t is successful 
iff 

hk{t)<xnixi{hk-i{t),hk+i{t)}. (6) 

Rule 2 {N = 2): At any site k where the update at- 
tempt was successful at (t — 1), at < we first randomly 
select a neighbor (left or right). This is equivalent to se- 
lecting either the left or the right border slice on the fcth 
processor. The update attempt is successful iff 

hk{t) < Kit), (7) 



where n is the randomly selected neighbor (n = k—1 for 
the left, n = fc -I- 1 for the right). At any site k where 
the update attempt was not successful at (f — 1), at < we 
keep the last n value. 

Rule 3 {N > 3): At any site k where the update attempt 
was successful at (i — 1), at t we first randomly select 
any of the N volume sites (indexed by rik) assigned to a 
processor. The selected site can be either from the border 
sites (either = 1 or Uk = N) or from the interior. The 
attempt is successful if the selected site is the interior 
site. When the border site is selected, the attempt is 
successful if condition {Tj) is satisfied. As in Rule 2, at 
any site k where the update attempt was not successful 
at (i — 1), a,t t we keep the last value. 

Rule 3 is essentially different from rules 1 and 2 in that 
it is a mixture of rule 2 and RD. At each t, depending on 
the selected Uk , the local update (deposition) at site k ei- 
ther follows rule 2 that requires checking with a neighbor 
or follows RD that just simply deposits a random positive 
real number rjkit). The probability (l—p) that the rule 3 
takes the form of RD is parametrized by N. At each t, the 
complementary probability p that rule 3 takes the form 
of rule 2 can be obtained by a direct count of lattice sites 
that have the assigned value either = 1 or = A'', 
and subsequently taking the configurational mean of this 
count. The mean density of these sites can be interpreted 
as the probability p{N) that at t a randomly selected site 
followed rule 2. In steady-state simulations (defined in 
Sec. imp this density does not depend on t and is found to 
be approximately p = ^JljN . This is because in simulat- 
ing conservative PDES the random selection of volume 
site rifc is not performed at every t for all lattice sites k. 
If at some t the update condition is not satisfied, the 
selection of n-^ must be postponed until some later time 
t' when condition ((TJ is satisfied. Explicitly, in the spirit 
of conservative PDES, if relation Q does not hold then 
the old rifc is kept for as many update steps as required 
until it finally holds at later t' . Note, if the draw of rik 
were performed at each t for all k the probability of se- 
lecting a border site would have been 2/A^, i.e., smaller 
than p. 

We define the utilization (u(t)) as the configurational 
average of the fraction of sites that performed an update 
at t. When A^ = 1, {u(t)) is simply the mean density of 
local minima of the interface. When N > 2, {u{t)) is the 
mean density of update sites. The velocity v{t) of the 
interface is defined as {dh{t)/dt). 

One distinction between the deposition models studied 
here and other restricted solid-on-solid models is that in 
the former the deposited random height increment 77 is 
a positive real number that can take on any value from 
an assigned real interval, while the latter usually con- 
sider integer height increments. Although in the context 
of applications to PDES this is the Poisson distribution 
with mean iip — I and unit variance that models the 
waiting times, represented here by 77, we also consider 
two other alternative depositions. One of them is uni- 
form deposition, where rj is sampled from the uniform 
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distribution with mean fijj = 1/2 and variance 1/12, re- 
stricted to the interval (0; 1]. The other one is a Gaussian 
deposition, where 77 is sampled from the Gaussian of unit 
variance and zero mean, restricted to the positive semi- 
axis. For this Gaussian, after nolmalization, the mean is 
fiG = ■\/2/7r and the variance is (1 — 2/7r). We find that 
varying the deposition type does not change the univer- 
sality class of the model. The purpose of introducing this 
variation is to better identify nonuniversal features in the 
initial evolution of the VTH interface. 



III. FLAT-SUBSTRATE CONDITION 

Rule 1 {N — 1), is a realization of random deposi- 
tion at local surface minima. As demonstrated in sim- 
ulations with Poisson depositions (56j . these interfaces 
belong to the KPZ universality class when the system 
size L is sufficiently large. Using standard finite-size 
scaling techniques Q for L of the order lO'^-lO^, the 
scaling exponents were determined numerically to be 
P = 0.326 ± 0.005 and a = 0.49 ± 0.01 |5|. A pecu- 
liarity of this scaling, present for all L, is the existence 
of the initial phase that does not scale, as illustrated in 
Fig-ffij. Another feature is the lack of scaling for small L, 
approximately L < 100. In this section we analyze the 
interface evolution for processes that obey rule 1 with 
Poissonian, uniform and Gaussian depositions, and show 
that the above lack of scaling is an artifact of the flat- 
substrate initial condition. 

Figure presents typical evolutions of the interface 



width for moderate to large L for both Poissonian and 
uniform depositions. A similar behavior is also seen for 
Gaussian depositions. A common feature is the existence 
of an initial growth phase, < t < to, where the widths 
do not scale. For t > to the widths obey the FV scaling, 
Eq. H2I3|) . with KPZ exponents. Figure [3 presents the 
FV scaling for Poissonian depositions (with 2a = 0.88) 
and for uniform depositions (with 2a = 0.94), both with 
a + z — 2, when the scaling transformation is applied for 
all t > 0. A similar picture of the data collapse is ob- 
tained for Gaussian depositions (with 2a = 0.92). The 
whisker-like structures in the growth part, clearly ob- 
served in Fig. 12 demonstrate the absence of scaling for 
< t < to- They vanish when the scaling is restricted to 
times t > to (the inserts to Fig. ^ and full data collapse 
is achieved for these later times. This initial transition 
phase is not a finite-size effect since to does not depend 
on L. For all L there is one common to that depends only 
on the deposition type. The largest to is observed for the 
deposition with the largest variance of the random height 
increments ry. In our examples, the smallest to is for the 
uniform depositions (variance 1/12) and the largest to is 
for Poissonian depositions (variance 1). In Gaussian de- 
positions (variance (1 — 2/7r)) the initial to falls between 
these two values. Thus, while the scaling shows that the 
mechanism of generating correlations (i.e., rule 1) has 
KPZ dynamics, the length to of this initial relaxation 
period to KPZ scaling is not universal. 

The initial transition period to is an artifact of the fiat- 
substrate initial condition. To investigate it further, we 
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FIG. 1: Time evolution of the simulated VTH width when 
V = 1 for various sizes L. (a) Poisson deposition, (b) Uniform 
deposition. The initial phase for t < to does not scale. 
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FIG. 2: Scaling of the VTH width when iV = 1. (a) Poisson 
deposition, (b) Uniform deposition. The inserts show the 
data collapse for t > to- The slope 2/3 of the growth phase is 
consistent with a + z = 2 and a = l3z. 
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FIG. 3: Time evolutions for A'' = 1 with Poisson deposition 
(L = 1000). (a) The widths, = f<w\ + />wS.. (b) The 
interface velocity v{t) and characteristic densities: the utiliza- 
tion {u(t)) {{u{0)) = 1, not shown), the simplex coefficients 
{f>{t)) and {f<{t)). The time to marks the transition to the 
steady state (the KPZ growth), and fx is the cross-over time 
to saturation. Ifere, fip — 1. 
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FIG. 4: Time evolutions for = 1 with Gaussian deposition 
{L = 1000). (a) The widths, w'^ = f<w% + ,f>w>. (b) The 
interface velocity v{t) and characteristic densities in analogy 
with Fig. 01 Here, fia = VV^- 



FIG. 5: Time evolutions for = 1 with uniform deposi- 
tion (L = 1000). (a) The widths, (b) The interface velocity 
v{t) and characteristic densities in analogy with Fig.^l Here, 
flu = 1/2. 



write out w'^{t) in its simplex form IgC 
linear combination 



as the convex 



(8) 



where f<{t) + />(t) = 1 is the convex sum, i.e., < 
f<{t)jf>{t) ^ 1- The characteristic densities f<{t) and 
/>(t) are the fractions of sites that have their heights 
less-then-or-equal-to and larger-than, respectively, the 
mean height h{t). The corresponding widths, computed 
on subsets that consist of these sites alone, are w%{t) 
and respectively. In individual simulations Eq.© 

is strictly satisfied and it is valid when averaged over 
many independent simulations. The convex sum is also 
valid for configurational averages of characteristic densi- 
ties. However, Eq.|(Sl) does not need to hold when char- 
acteristic densities and widths are changed to their cor- 
responding configurational averages (because, in general, 
(ab) 7^ (a)(5)). Configurational averages of characteristic 
widths and densities, and the interface velocity v{t) are 
presented in Fig. Ofor Poissonian depositions, in Fig. 0| 
for Gaussian depositions, and in Fig.[51for uniform depo- 
sitions. At t = 0, the interface velocity and the utilization 
have their highest values (u(0)) — 1 because Eq. © is 
satisfied at all sites. This first step at t = is simply 
a random deposition step. The mean height h(0) — fi 
is the mean of the distribution from which rj is sampled, 
which is lip — 1, fiG — \/2/7r and fj,u = 1/2 for Pois- 
son, Gaussian and uniform depositions, respectively. The 
fraction /> (0) of sites that have their heights larger than 
h{0) is easily computed from the corresponding distribu- 
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tions as the probability of selecting a site that has h{0) 
larger than fi. This gives for Poisson deposition />(0) — 
J^^ dx exp{—x) — 1/e 0.367; for Gaussian deposition 

MO) = y^QdxeM~^V2) - l-crf(l/V^) « 

0.428; and, for uniform deposition />(0) = /^^^ dx — 1/2. 
These fractions and their complements /< (0) = 1 — /> (0) 
are clearly observed in Figs. I3I5I Correlations between 
lattice sites start to build up at t = 1. Since initially the 
density f<{t) is larger than f>{t), depositions take place 
more often at sites with h <h then at sites with h > h. 
This causes f<{t) to fall and />(i) to rise (Figs. |3>©) 
and a faster growth of 'w%{t) than Wy{t) (Figs. |2h-Eh)- 
On the average, as the density (/> (t)) rises the density of 
local minima {u(t)) decreases. This initial evolution from 
the RD surface at t = to a surface with correlations at 
to ends when (/>(i)) « (/<(*)) « 1/2. As Figs. EE il- 
lustrate, at to the simulations attain what we label as a 
steady state, one that is characterized by a constant uti- 
lization. We show in Sec.^that v{t) is related to {u{t)) 
by a simple linear relation, hence, the steady state can 
be alternatively defined by a constant velocity. 

The correlated growth phase when the scaling is ob- 
served, when to < t <^ tx, is characterized by a slight 
but noticeable excess of (/>(t)) over {f<{t)). At satu- 
ration, for t ^ tx, (/>) ~ (/<)■ The densities {f<{t)) 
and (/>(t)), and the widths (w<(t)) and (w>(t)), provide 
the information about the height distribution ^{h/h) of 
the interface local heights about the mean height h{t). 
It is transparent from Figs. I3I5I that for early times, 
t < to, ^{h/h) is characterized by a positive skewness 
and evolves to approximately a symmetric distribution 
at to- This distribution function for Poisson depositions 
at local minima is presented in Fig. |S1 The computation 
of $(a;) is outlined in the Appendix. 
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FIG. 6; Poisson deposition at local surface minima for L — 
1000: distribution function ^{x) of the interface local heights 
X — h(t)/h{t). (a) Early times t < to- (b) The steady state 
for t > to. Cubic-spline curves through the simulation data 
(symbols) are guides for the eyes. Here, to oc 100. 



The skewness of the height distribution in the station- 
ary state of the KPZ growth has been analyzed before by 
den Nijs and co-workers |U 163. l65l | for Kim-Kosterlitz 
models with integer step-height differences. They report 
that KPZ scaling is realized at times larger than a char- 
acteristic time scale that is related to slope densities. In 
(1-1-1) dimension the KPZ dynamics is characterized by 
zero skewness because the height distribution of the sta- 
tionary state is Gaussian [6^ . In our models, the growth 
can be characterized alternatively either by the density 
of local minima (i.e., the utilization for = 1) that is 
the same as the density of local maxima [s^] or by the 
density of local slopes. Since all these densities sum up to 
one, the constant utilization in our model is equivalent to 
a constant density of local slopes. Explicitly, we define 
the steady evolution state (or the steady state simula- 
tions) as the evolution that has the following character- 
istics: 1. the density of update sites is constant; and, 2. 
the skewness of the height distribution is approximately 
zero. Starting from the flat substrate, the steady growth 
state is achieved after the initial relaxation time to- In 
the steady state the KPZ scaling is clearly observed. 

The initial time interval from t = 1 to t = to can be 
interpreted as the time scale over which the system re- 
tains the memory of the flat-interface initial condition. 
This time is a nonuniversal parameter that depends on 
the variance of the distribution from which the random 
height increment rj is sampled. The existence of this time 
scale accounts for the absence of universal scaling for 
small system sizes, even if the rule that simulates the 
growth represents a generic KPZ process. For KPZ dy- 
namics the characteristic time scale on which the corre- 
lations are being built is of the order of the system size 
tx ~ L^/^. If this time scale is smaller than the memory 
scale, L^/^ < to, the interface saturates before the simu- 
lations reach the steady state, and for such L the KPZ 
scaling is not observed. The universal KPZ scaling is 
clearly observed when the system size is large enough to 
loose the memory of the initial condition on time scales 
smaller than L^^^, i.e., when to ^ ■ 



IV. SCALING ANALYSIS 

In this section we analyze the interfaces generated by 
deposition rule 3 that represents two simultaneously act- 
ing growth mechanisms: one is RD and the other is de- 
position rule 2, both with Poisson-random height incre- 
ments. First we show that rule 2 generates KPZ correla- 
tions. Then we obtain the universal scaling function for 
interfaces produced by rule 3. 

Although rule 2 allows the /cth site to accept a deposi- 
tion even if it is not a local minimum, this rule has all the 
essential characteristics of rule 1, examined in Sec. lIIII At 
each time step a site must compare its local height with 
a local height of at least one of its immediate neighbors. 
As in rule 1, deposition may not happen at a local max- 
imum. But, since now it may happen either at a local 
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FIG. 7: Time evolutions for Poisson deposition when N = 2 
and L = 1000. (a) The widths, (b) Characteristic densities 
and the interface velocity, in analogy with Fig. |3 
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FIG. 8: Scaling for t > to when N = 2. The initial relaxation 
period for t < to is clearly noticeable in the insert that shows 
the widths before scaling. 



minimum or at a local slope, the utilization of rule 2 is 
larger than the one of rule 1 (compare Fig.|3) and Fig.|3)) 
so the interface velocity is larger. Other than that there is 
no difference between these two deposition mechanisms, 
and the analysis presented in Sec. Illll for interfaces pro- 
duced by rule 1 can be restated for the interfaces that 
grow by rule 2. In particular, both growths evolve on the 
same time scales (Fig. [T)), with the initial memory scale 
to oc 100. Figure |H1 shows the scaling function for the 
interface widths for t > to, obtained with 2a = 0.9 and 
z = 2 — a, characteristic for KPZ scaling. Thus, these 
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FIG. 9: Time evolutions for Poisson deposition when A'' = 
100 and L = 1000. (a) The widths, (b) Characteristic den- 
sities and the interface velocity, in analogy with Fig. |3] and 
Fig. El 



interfaces belong to the KPZ universality class and the 
scaling function is given by Eqs. (|2I3|I . A small departure 
of 2a from one, also present for N ~ 1, is discussed in 
Sec lyTBl 

Deposition rule 3 produces a larger utilization than 
rule 2 because now, depending on TV, deposition at site 
k may sometimes be accepted regardless of the relation 
between its local height and local heights of its neigh- 
bors. Now at each t, any site may increase its height, 
including a local maximum. Probability p{N) that a 
site has to compare its local height with a neighbor is 
the probability of applying the rule 2, which is the only 
mechanism that creates correlations. Alternatively, de- 
position at the k site may happen as RD. A combination 
of these two deposition mechanisms produces a similar 
time evolution of characteristic densities and the widths 
as it is observed when rule 2 is acting alone, except that 
now the transition to steady-state simulations and the 
cross-over to saturation take place on larger time scales 
(compare Fig. [71 and Fig. O. In particular, the initial 
lack of scaling extends to to{N) cx toN/2, where marks 
the end of the initial relaxation period in the worst-case 
scenario simulations. This initial relaxation time to{N), 
when the system "remembers" the flat-interface initial 
condition, manifests itself in the evolution of interface 
widths as an early growth phase fFig. llOl) that follows the 
RD power law with /3i = 1/2. The later growth phase, 
to{N) <^t<^tx, follows the power law with /32 = 1/3 — e, 
where e is a small positive number. The evolution of the 
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FIG. 10: Two growth phases in the time evolution ol the 
width, simulated with rule 3 for N > 3 and L = 1000. 



interface width can be summarized as 




t < to{N) 

to(A^)<<<ix(iV) 
t>ix(iV), 



(9) 



where g{N) is a monotonic function of TV. After perform- 
ing scaling in L of the saturated width, it appears that 
g{N) is Unear. The first growth phase is the initial relax- 
ation when {w'^{t)) does not scale, therefore the following 
analysis is valid only for t > to{N). 

From the point of view of scaling, {w'^{t)) is a family of 
curves parametrized by L and N. Fiaure lTTt t presents the 
saturated width {w^) plotted against N for selected val- 
ues of L. These curves can be scaled in L so as they col- 
lapse onto one curve. Figure ITTt) shows the scaled width 
where 2a = 0.9. Since Iogio((w2)/L2a) _ 
logiQ N + const, (dashed hne in Fig. \llb). values at sat- 
uration may be further scaled in N, which gives the col- 
lapse to one point {w^) / {N L"^^) k, const. The order of 
scaling can be reversed, i.e., scaling in N can be followed 
by scaling in L, leading to the same result. Globally, the 
saturated width plotted vs {NL) follows a line of slope 
one fFig. ITTb'). 

Time evolution of the scaled width {w"^ (t)) / {N L^") for 
t > to{N) is displayed in Fig.^l Here, to obtain a per- 
fect align at saturation with the curves for TV = 2, which 
facilitates further scaling in t, we introduced a small cor- 
rection such that NL'^" is multiplied by (lie), where 
e is a small fraction (explicitly, e is the relative spread 
of data about the fit at saturation; the scaling is clearly 
seen with e = 0). Since the only mechanism that induces 
correlations in this model is deposition rule 2, which pro- 
duces surfaces from the KPZ universality class, it is ex- 
pected that for each N the length scales should couple 
with time scales via a dynamic exponent that satisfies the 
KPZ identity z — 2 — a. Indeed, as the partial scaling 
in t with respect to L shows (Fig. \l'6b). this choice gives 
a good data collapse for each N and this is not changed 
noticeably when z = 1.5 is chosen. In Fig. 113b . read- 
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FIG. 11: The saturated interface width as a two-parameter 
family of curves obtained in simulations with rule 3: (a) as 
a function of A'^ volume elements per lattice site for several 
lattice sizes L; (b) as a function of A*', after scaling in L of 
the data in figure (a) (the dashed line of slope one is plotted 
as a reference); (c) data from figure (a) plotted vs (NV) to 
see global trends. The data align along a straight line of the 
mean slope one (the dashed-line envelope). 
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FIG. 12: The steady-state evolution for t > to of the scaled 
widths obtained with rule 3 as a two-parameter family of 
curves for: (a) N = 10; (b) iV = 100; (c) L = 100; and, 
(d) L = 1000. 



ing from the left, the groups represent iV = 2, = 10, 
iV = 100, and iV = 1000. The last step is the scaling in 
t with respect to N of the results displayed in Fig. W3a . 
Here the inspection of the dilatation with N of the ini- 
tial ^0 proves useful since it leads to the observation that 
for any h the initial to{N) can be expressed approxi- 
mately as {N/2)to{N = 1 or 2). The transformation 
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FIG. 13: Scaling for the curves of Fig. [H (a) in t with 
respect to L, with the KPZ dynamic exponent z = 2 — a.\ (b) 
in t with respect to A'^ of the data in figure (a). Data labels 
displayed in the legend are common for both figures (a) and 
(b). 



t t/{N/2) shifts (to the left) all the curves for > 2 
into one position. The family of curves for A'^ = 2 is 
shifted into this position when t t/2. The final result 
is the scaling function shown in Fig. 113b. The mean slope 
of the growth part is 2(3 k, 0.58 ± 0.02, consistent with 
the slope obtained for 2a = 0.9 from the KPZ relation 
2/3= (2a)/(2-a) «0.58. 

The final result, valid for t > to{N), can be summa- 
rized as 



(10) 



where f{y) satisfies Eq. (jJJ, and z = 2 — a with a « 
1/2. Accordingly, the interfaces generated by the deposi- 
tion/update rule 3 belong to the KPZ universality class. 
In the scaling regime, the evolution can be written out 
explicitly as 



NT2a 



i >tx(iV), 



where ty,{N) - {N/2)L^ 



(11) 



V. INTERFACE VELOCITY 

For the deposition/update models considered in this 
work it is possible to find the exact relation between 



the utilization and the interface velocity. The velocity 
v{t) is defined as the configurational average of v{t) = 
dh{t)/dt — {^/L)J2k=id-hk{t)/dt. Translating the up- 
date recipe to a continuum version with a continuum 
time-step increment St, the update operation at site k 
can be summarized as 

hk {t) + St rjk (t) , on update 
/ifc(i) , otherwise. 



hk{t + St) 



(12) 



Substituting the above to the definition of dhi;{t)/dt as 
the limit when 5i — > 0, gives at t 



dhkjt) 
dt 



rjk{t) , k is the update site 
, otherwise. 



(13) 



In the set of L sites the number of update sites is M{t) = 
Lu{t). Since only at these sites dhk{t)/dt is not zero, the 
mean is 

M(t) M(t) 



1 E 



u{t) 



1 



k=l 



M{t) 



(14) 



fc=i 



Let /ip be the mean of the distribution P{vi) from which 

rjk is sampled. In the limit of M ^ oo, (1/Af) X]/c=i '^fc ~^ 
J^drif] P{r]) = Hp. Thus, for sufficiently large M, the 
second factor in Eq. p4|) is /ip = const. Taking the 
configurational average of Eq. 114|l gives 



v{t) = {u{t)) Hp. 



(15) 



The above derivation can be repeated for the discrete 
case, taking St — I. Equation (|15|l is strictly satisfied by 
the simulation data fFigs. 171 and . In simulations 
v{t) is computed numerically with St = 1, v{t) = h{t) — 
h{t— 1), and v{t) is obtained by averaging v{t) over many 
independent simulations. 



VI. DISCUSSION 

In our steady-state simulations for A^ > 1, the mean 
density of sites p{N) = yj2/N can be interpreted as the 
probability that at t a randomly selected site followed 
rule 2. The complementary density p = 1 — p is the 
probability that a randomly selected site followed RD. A 
similar interpretation can be given to the mean density 
of update sites {u(t)) as the probability that at t a ran- 
domly selected lattice site increased its height; however, 
{u{t)) does not define a probability distribution [HtI IH^. 
Using a recently introduced discrete-event analytic tech- 
nique for surface growth problems [H^, El] , it is possible 
to derive a mean-field-like expression for the utilization 
(m) in the steady-state simulations. For N > 2 and L > 3 
this expression ;58j is 



= 1- 



p{N) 



1 - 



p{N) L - 1 



(16) 



Because {u{t)) and v{t) are related by a constant mul- 
tiplicative factor /j,p, by Eq. (|15|l . in the scaling regime 
for t > to{N), Eq. (|16|) also gives the interface velocity. 
However, for early times {u(t)) is unknown. 
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A. Kinetic roughening 

The scahng expressed by Eq. H10() can be written in a 
more general form by incorporating the probability p — 
y/2/N. This gives 



L^WfpH 



(17) 



The analysis of Sec. llVl can be repeated for {w{t)), which 
gives equivalently 



(18) 



The transition time to the scaling regime is io(p) = ^o/p^ ^ 
where tg marks the end of the initial non-scaling period in 
the case of simulations with the deposition/update rule 
2 acting alone. The cross-over time to saturation ty, can 
be read directly from Eqs. p7ll8|l . p^ty, {p)/L'' « 1. This 
gives tx{p) ~ t-xjp^, where ix is the saturation time 
when rule 2 acts alone. 

Our results for the scaling function, Eqs. I|1(JI11|I and 
Eqs. (|17I18|) . are gener ally in accord with the work of 
Horowitz and Albano [Sg and Horowitz et al. [s^ . 
who analyzed (1-1-1) dimensional two-component solid- 
on-solid models mixed with RD. In Ref. 38], the growth 
is simulated by ballistic deposition (of the KPZ univer- 
sality class) that takes place with probability p and by 
RD that happens with probability {1 ~ p). In Ref. |33 
the deposition model mixes RD (taking place with prob- 
ability (1 — p)) and random deposition with surface re- 
laxation (of the EW universality class) that takes place 
with probability p. A common characteristic of these 
two models and our model is the presence of two growth 
phases in time evolution of the interface width (Fig. ^| 
and Eq. 0), where the early phase follows RD growth 
and this is the phase that defies the universal scaling (fig- 
ures 4 in Refs. H^ill^)- Although this initial absence of 
scaling is not studied in Refs. |37ll38l |. based on our anal- 
ysis of Sec. mil wc infer that this early RD growth seen 
in Refs. |37l |38j must be a long-time effect of some par- 
ticular initial condition (not stated in ^3^, Isll) adopted 
in these simulations. As we showed in Scc. lIIII the length 
of this initial memory scale is a nonuniversal parameter, 
so comparing cross-over times to the scaling regime does 
not contain useful information. But, the existence of this 
initial memory scale and the existence of scaling for times 
larger than the initial relaxation time are both universal. 
The scaling law reported in Refs. j^ll^ is 



W{t,L,p)^-^Fi^^ 



(19) 



where the exponents a and z are characteristic for the 
universality class of the model (i.e., KPZ in [s^l and EW 
in )■ For the mixture of RD and random deposition 
with surface relaxation (the EW type model) the numer- 
ics give 5 « 1 and y k, 2 j37j, in agreement with our 



findings, Eq. (|18|l . However, for the mixture of RD and 
ballistic deposition (the KPZ type model) it is conjec- 
tured in |33| that (5 « 1/2 and y « 1. In this latter 
scaling, although the relation between y and S, y — 2d, 
agrees with our findings, these powers are by the factor of 
2 smaller than our values. As we analyze later in this sec- 
tion, this difference implies that according to Ref. [s^ 
the admixture of RD should affect time scales as pt, while 
our findings clearly indicate the p^t behavior. 

The sensitivity of the surface evolution to the initial 
condition has been recently pointed out by Kortla et al. 
32] in relation with phase ordering in two-component 
solid-on-solid models. In our modeling, to incorporate 
fully the dynamics of the growth, the KPZ equation 
should contain the mean interface velocity v{t) 



ht = v{t) +iyh^^ + ^hl + C 



(20) 



Equation can be valid only for the steady-state sim- 
ulations, when v(t) = const. Therefore, it does not de- 
scribe the initial dynamics at t < to of our deposition 
models, while Eq. H20() does. Results obtained in Sec. lIVI 
summarized by Eqs. H17I18|I . indicate that the admixture 
of the RD process (present with probability (l—p) in rule 
3) elongates time scales in Eq. ipU]) . We now scrutinize 
this effect from the point of view of the affine transfor- 
mations involved. 

Consider the following transformations that are ex- 
pressed by Eqs. H17I18|I 

x^x' h^h' = ^h, t^t' = ^t, (21) 

where g{p) is to be determined (note, g{p) is not inde- 
pendent because x, h and t are transformed simultane- 
ously). Denoting v — —h^, in (1-1-1) dimensions the con- 
vective derivative in Burger's flow is Dtv — Vt + Xvvx |3] ■ 
It is straightforward to show that DtV is invariant (up 
to a multiplicative factor) under transformations (|21|l if 
a + z = 2, g{p) = p and if A ^ A', Ay = A. Then w^{t) 
transforms to w''^{t'), w''^{t') = w'^ {t) / {L'^" / p'^) ~ $(t') 
or, equivalently, w{t) w{t)/{L"/p) ^ ^{t')- This gives 
Eqs. P7I18|I . If we set A = in Eq. ^ ( assuming the 
scaling regime of constant v{t)), then the resulting EW 
equation is invariant under the same transformations pro- 
viding z = 2, 2a + 1 = z and v v' , v'p^ = v. This, 
again, gives Eqs. H17I18|I but with z — 2 and a = 1/2. 
Transformations defined by Eq. (|21|l can be seen as su- 
perpositions, where the first scaling 

x^x/L, h-^h/L", t^t/L^ (22) 

is followed by 

X— s-x, h ^ ph , > p^t . (23) 

(The order in which transformations H22|l and H23|) are su- 
perimposed to form transformation (|21|) can be reversed.) 
The second transformation, defined by Eq. H23|) . leaves 



12 



the KPZ equation (|20|l invariant provided that new co- 
efHcients v' and A' are related to the old coefScients by 
the relations v']? — v and A'p^ = A (which is the result 
obtained before) and the old velocity v(t) changes to a 
new velocity = v{t')/p. Notice, in the absence of 
the RD admixture we have p — 1 and the transforma- 
tion (|23|) is the identity. Then transformation H21|) is the 
usual FV scaling, Eqs. (|2I3|I . with proper choices for a 
and z, depending on the process at hand (either KPZ or 
EW). When the RD process is present we have p < 1. 
Then, inverting transformation (|23|l gives the changes in 
the local-height field h/p and in time scales t/p^. Since 
the elongation of the time scale is inversely proportional 
to p^ this effect is more pronounced than the stretching in 
h. These phenomena take place for alH > since trans- 
formation H23|l is valid for all i, while the FV scaling or, 
equivalently, affine mapping given by Eq. takes place 
only in steady-state simulations at times t 3> to(?') when 
the system has lost the memory of the initial condition. 

One possibility for the behavior seen in this paper for 
an admixture of the KPZ fixed point with the RD fixed 
point could be for a reason similar to the floating-fixed 
point seen in critical phenomena [6^ IbtI l6^ . In that 
case two physical fixed points (in integer dimensions) are 
joined by a line of fixed points that are inaccessible (in 
that case in non-integer dimensions). The result is that 
for a finite system size, and depending on the boundary 
conditions involved, there is an effective critical exponent 
[HlzSEllli- Only in the limit of infinitely large sys- 
tems is the physical fixed point approached. The effective 
critical exponents along this line of fixed points satisfy 
(hyper)scaling relationships, just as here the relationship 
a -\- z — 2 holds. Further investigations studying much 
larger systems would be needed to see if the floating- 
fixed point picture holds in our case where properties of 
both the RD and KPZ fixed points are mixed into the 
nonequilibrium surface model. 

In summary, the only effect of the RD admixture to 
either a genuine KPZ or EW process is the dilatation of 
growth scales. The RD blending does not change the uni- 
versality class of the interface since it does not change the 
dynamics of mechanisms that are responsible for building 
correlations. However, such dilatation, when combined 
with the initial flat-substrate condition, may obscure a 
clear observation of KPZ scaling in simulations as well 
as in experiments. 



B. Roughness 

In this section we discuss flnite-size effects in scaling for 
the VTH interfaces simulated by the worst-case-scenario 
deposition/update rules 1 and 2. Despite that these in- 
terfaces belong to the KPZ universality class, the precise 
value of the roughness exponent depends on the lattice 
size L and on the type of deposition. In the case of Pois- 
son deposition it requires large L to attain scaling with 
2a « 1 (a approaches 1/2 from below as L is in- 




FIG. 14: Distributions $(x) of the interface widths at sat- 
uration, X — /{w^). Results of simulations (symbols) are 
compared to Eq. 1)24^ (continuous curve), (a) Uniform de- 
position when L = 1000 and A'^ = 1,2 (binning interval 
Au = 0.05; the data for A'^ = 1 fall on the top of data for 
N — 2). (b) Uniform and Poisson depositions when L — 1000 
and iV = 1,2 {Au = 0.05, Ap = 0.5). (c) Poisson deposition 
when L = 1000, 10000 and Af = 1, 10, 100 (Ap = 1). 



creased), while for uniform depositions the required L is 
approximately two orders smaller than that for Poisson 
depositions. 

At saturation the growth of a (l-fl) dimensional KPZ 
(or EW) surface can be mapped onto a diffusion problem 
with column- height fluctuations Sh ~ ^/z^, where z^ 
is the dynamic exponent of a random walker 0, Q that 
connects with the roughness exponent, az^ = 1. The 
exact a = 1/2 indicates the total lack of correlations 
and a > 1/2 indicates their presence. The random- walk 
interfaces are characterized by the following width distri- 
bution function ^(w^ j (w^)) 18] 



2 ^ 



n exp [ — g""' 



(24) 



In simulations of growth models where depends on a 
single length scale, $ is obtained by normalizin g a his- 
^ram PiyJ^) of the width distribution HI [zl H El 



= (W2)p(u;2). 



(25) 



This technique of describing surfaces in terms of their 
scaling functions $ gives very good agreement between 



13 



theoretical functions <I> (whenever available) and simu- 
lated curves for growth models with integer height in- 
crements [H m [H m [l^. Figure [H shows $ ob- 
tained in simulations with several N and L for our depo- 
sition/update models with Poisson and uniform deposi- 
tions. We observe that these curves closely follow the the- 
oretical curve given by Eq. H24|) . In the computation of 
the quantities in Eq. (|25|l (see Appendix) we used a vari- 
able step size A in binning the data (A = 0.05, 0.1, 0.5 
and 1) and a variable number N^ata of data points at 
saturation (6.4 x 10^ < N^ata < 10*) to ensure that the 
results of Fig. E| represent the true limit of <& obtained 
in our models. 

The exact collapse of the distributions <& obtained in 
simulations on the theoretical curve of Fotlin et al. |0| , 
given by Eq. 1)241) . is not related to the type of the de- 
position, as is illustrated by the data in Fig. E| and is 
observed even for small system sizes (L oc 100). There- 
fore explaination for the sensitivity of a to L and to the 
deposition type must lie elsewhere. For moderate to large 
L presented in this work, when = 1 data collapse at 
saturation required 2a = 0.94 for uniform deposition, 
2a = 0.92 for Gaussian deposition, and 2a ^ 0.88 for 
Poisson deposition. Similarly, when N — 2 the collapse 
was achieved with 2a = 0.94 for uniform deposition and 
2a = 0.9 for Poisson deposition. This variation of a 
with the deposition type suggest that a small departure 
Sa = l — 2a from the exact value 2q! = 1 is a nonuniversal 
parameter. As observed, the largest difference Sa is for 
depositions that have the largest variance (i.e., Poisson) 
and the smallest Sa is for depositions of the smallest vari- 
ance (i.e., uniform). This observation strongly indicates 
that Sa is related to the system memory (introduced in 
Sec. i.e., to time scales T on which the interface 

does not remember past depositions. In other words, Sa 
depends on the minimal interval T such that, on the av- 
erage, a deposition at time t does not affect depositions 
at time {t + T). Long memory scales T lead to a build- 
up of temporal correlations, thus, produce a variation 
to Gaussian noise in Eq. possibly modifying noise 
strength D in Eq. Such a variation may influence 
the way in which the system attains the stationary state 
and, possibly, the value of a. We leave this issue 
open to future investigations. 
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FIG. 15: The false-scaling time evolution for t > of the 
scaled widths, obtained with Poisson depositions in rule 3: 
(a) iV = 100; (b) N = 1000; (c) L = 1000; (d) L = 1000. 
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C. False scaling 

A blind application of the scaling technique to numer- 
ical data sets, without supporting analysis of temporal 
scales and invariants involved, may lead to false conclu- 
sions. For example, assuming a global linear behavior of 
(w^) vs the variable {NL) fFiEclllfc). it is possible to col- 
lapse the saturated widths by using an effective exponent 
2a « 1. Performed for all t > 0, such scaling shows two 
distinct growth regimes in the evolution curves (Fig. I15|l . 
For each L, these can be further collapsed into groups by 
scaling t t/N""^ for aU t, taking zi = 1 (Fis.HRi). The 



FIG. 16: False scaling for the curves of Fig. 1151 (a) Scaling 
with A'^ in t, zi = 1. (b) Scaling with L in t for the data in 
figure (a). 



excellent total data collapse in the first growth regime (of 
slope 2/3i in Figs. ll5ll6jl is obtained when t -> t/{NLY\ 
but at the expense of only an approximate collapse in the 
second growth regime (of slope 2/32 in Figs. 11511^1) . This 
final result (Fig. ITBt )') seems to look plausible because 
2[3i — 1 and 2a = 1 give a hypothetical exponent zi = 1, 
and the curves fall exactly on the top of each other in the 
early growth phase and at saturation. However, upon a 
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closer inspection it appears that the length of the inter- 
mediate regime of slope 2(32 expands when L is increased. 
In the limit of large but finite L, the scaled curves 
in Fig. 116b will cover the lower right-hand semi-plane, 
bounded by lines y = \ogi^{{'uP') / {N L)'^°') — const, and 
y — 2/3i log]^Q(t/(Afi)^i) + const. Hence, there is no scal- 
ing. In fact, this example supports our earlier conclusion 
that the initial growth phase for t < to{N), the artifact 
of the initial condition, does not scale. 



D. Application to PDES 

In modeling of the conservative update mode in PDES, 
we represent sequential events on a processor in terms 
of their corresponding local virtual times. The column 
height that rises at the kth lattice site in the simulated 
VTH represents the total time of operations performed 
by the fcth processor. These operations can be seen as 
a sequence of update cycles, where each cycle has two 
phases. The first phase is the processing of the assigned 
set of discrete events (e.g., spin flipping on the assigned 
sublattice). This phase is followed by a messaging phase 
that closes the cycle, when a processor broadcasts its 
findings to other processors. But the messages broad- 
casted by other processors may arrive any time during 
the cycle. Processing related to these messages (e.g., 
memory allocations/deallocations, sorting and/or other 
related operations) are handled by other algorithms that 
carry their own virtual times. In fact, in actual sim- 
ulations, this messaging phase may take an enormous 
amount of time, depending on the hardware configura- 
tion and the message processing algorithms. In our mod- 
eling the time extent of the messaging phase is ignored 
as though communications among processors were taking 
place instantaneously. In this sense we model an ideal 
system of processors. The local virtual time of a cycle 
represents only the time that logical processes require to 
complete the first phase of a cycle. Therefore, the spread 
in local virtual times represents only the desynchronisa- 
tion that arises due to the conservative algorithm alone. 

The measure of this desynchronization is provided by 
{w{t)). Since in PDES the memory request per processor, 
required for past-state savings, is determined by the ex- 
tent to which processors get desynchronized, {w{t)) can 
be considered as an indirect approximate measure of this 
memory request. Its growth during the entire time span 
of the PDES computations is given explicitly by Eq. ^ . 

We showed that in conservative PDES, given the PDES 
size (finite L and A^), this memory request does not grow 
without limit but varies as the computations evolve. The 
fastest growth, proportional to \/t, characterizes the ini- 
tial start-up phase. The length of the start-up phase de- 
pends on the load per processor (represented here by N) . 
The start-up phase is characterized by decreasing values 
of both the utilization and the progress of the global sim- 
ulated time (i.e., the smallest local virtual time from all 
processors at each simulation step). In the steady-state 



simulations, when the utilization has already stabilized 
at a mean constant value (and so has the progress of 
the global simulated time), the memory request grows 
slower, at a decreasing rate ~ In this phase, the 

mean request can be estimated globally from Eq. H10(l 
or Eq. Hll() . The important consequence of scaling, ex- 
pressed by Eq. Hl()|l . is the existence of the upper bound 
for the memory request for any finite number L of pro- 
cessors and for any finite load N per processor. As it 
is stated by Eq. Hll|l. on the average, this upper bound 
increases proportionally to \/NL with the size of conser- 
vative PDES. 

The characteristic time scale to{N) from the first step 
to the steady-state simulations can be estimated by mon- 
itoring the utilization for the minimal processor load (to 
determine to) and, subsequently, scaling this time with 
N. Similarly, the characteristic time scale to tx{N), 
when the desynchronization reaches its steady state, can 
be scaled with the processor load to determine an ap- 
proximate number of simulation steps to the point when 
the mean memory request does not grow anymore. 

During the steady state simulations the utilization (u) 
is given by the approximate Eq. H16|) . where p{N) — 
y'2/N. The smallest utilization {up) for any processor 



load is obtained by taking the limit L 



of Eq. (Uni: 



{up) = {l-p{N)/2){l-p{N)/4). 



(26) 



Since as TV — > cx), p{N) decreases as ^ Eq. 
shows that (up) grows very fast when the processors' load 
increases. For example, (up) for N = 100 is about 90%, 
close to its asymptotic limit of 100%. For the minimal 
processors' load {p — 1) Eq. ((TB|) gives for L > 3 



3 1 



(27) 



In the limit L — > oo, (uq) is the smallest possible value 
of the utilization. As Eq. (|27|l shows, this value is a 
non-zero constant (equal to 3/8 when derived from sim- 
ulations with Poissonian distribution of waiting times). 
This nonzero lower bound on the utilization and the fi- 
nite upper bound for the memory request for finite L 
show that conservative PDES are generally scalable with 
the number of computing processors when performed in 
the ring communication topology. The extension of this 
conclusion to other communication topologies requires a 
separate study. 



VII. SUMMARY AND CONCLUSIONS 

We considered a two-component growth in (1+1) di- 
mensions. One of the components is RD that takes 
place with probability (1 ~ p). The other component, 
which takes place with probability p, is a deposition pro- 
cess that generates correlations typical of KPZ dynamics. 
The growth is simulated from an initially flat substrate. 



15 



We show that the flat-substrate initial condition is 
responsible for the existence of the initial non-scaling 
regime in simulations. The length of this initial phase 
is a nonuniversal parameter (it depends on the type of 
depositions and on the particulars of the model). How- 
ever, its presence is a universal phenomenon. 

During the initial phase the simulations relax to a 
steady state. For the models considered in this work, 
the transition time to the steady state can be defined as 
the time when the mean interface velocity attains a con- 
stant value. We showed that for these models the mean 
interface velocity is a multiple of the mean utilization. 

During the steady state the interface width satisfies 
FV scaling. We derived the universal scaling function for 
the width and showed that the RD admixture acts as a 
dilatation mechanism to the time and height scales, but 
leaves the KPZ correlations intact. This conclusion has 
been generalized to two-component models that mix RD 
with depositions that classify within the EW universality 
class. In particular, we showed that the RD admixture 
is responsible for the p dependent affine change of scales 
(h — > h/p and t t/p^) that is superimposed on the 
usual scaling and leaves the dynamics invariant. 

The models, studied in this work, that give rise to 
the KPZ correlations are the Poisson-, the Gauss-, and 
the uniform-random depositions either to local interface 
minima or to local minima and randomly selected lo- 
cal slopes. Despite that the simulated interfaces be- 
long to the KPZ universality class, the precise value 
of their roughness exponent depends on the deposition 
type. This observation suggests that such noisy deposi- 
tion mechanisms may produce relatively long-scale tem- 
poral correlations. Secondly, this small departure from 
the exact value is nonuniversal. Further studies are re- 
quired to investigate this issue. 

In application to conservative PDES, we showed that 
the memory request per processor, required for state sav- 
ings, does not grow without limit for a finite number of 
processors and a finite load per processor but varies as 
the PDES evolve. The important consequence of the de- 
rived scaling is the existence of the upper bound for the 
desynchronization, thus, for this memory request. Also, 



the utilization of the parallel processing environment has 
a non-zero lower bound as the number of processors in- 
creases infinitely. Thus, the conservative PDES are gen- 
erally scalable in the ring communication topology. 
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APPENDIX: DISTRIBUTIONS 

Both hk{t) and w'^{t) are real numbers that take on 
continuous values. Suppose S = {/i, /2, /at} is a 
set of N such numbers obtained in simulations. Let 
a — min(S') and b = max(S'). The interval {b— a) is par- 
titioned into M segments, each of length A = {b~a) /M. 
Each segment is a bin that is indexed by its left end 
Hi = a + {i — 1)A, 1 = 1,2, M. Let rrii be the mul- 
tiplicity of the zth bin, i.e., m.i is the number of points 
from S that fall between yi and The mean value 

of S is {f) = l/A^Ef=i fk = EZi P^ where (/.) = 
{l/rrii) Y^^Li fin is the mean value taken on a subset of S 
that belongs to the ith bin, and Pi is the frequency func- 
tion Pi = {mi/N){l/ A). The abscissas Xi and function 
values $j = $(a;,) are: Xi = {fi)/{f) and $i = {f)Pi. 
These values are plotted in Fig. and Fig. ^1 Here, 
^{x) is properly normalized dx^{x) = 1. The abso- 
lute spread (b — a) determines a suitable A in the compu- 
tation of (fi) within acceptable precision. This gives the 
total number M of bins. For a given M, the accuracy of 
Pi depends on N. 
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